Phase space and quark mass effects in neutrino emissions in a color superconductor 
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We study the phase space for neutrino emissions with massive quarks in direct Urea processes in 
normal and color superconducting quark matter. We derive in QCD and the NJL model the Fermi 
momentum reduction resulting from Fermi liquid properties which opens up the phase space for 
neutrino emissions. The relation between the Fermi momentum and chemical potential is found to 
be pf ~ /tt(l — k) with k depending on coupling constants. We find in the weak coupling regime 
that k is a monotonously increasing function of the chemical potential. This implies quenched 
phase space for neutrino emissions at low baryon densities. We calculate neutrino emissivities with 
massive quarks in a spin-one color superconductor. The quark mass corrections are found to be 
of the same order as the contributions in the massless case, which will bring sizable effects on the 
cooling behavior of compact stars. 
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Ch ' I. INTRODUCTION AND CONCLUSION 

3 ' 

The baryon density in the core of a compact star is likely to reach several times the nuclear saturation density, po ~ 

0. 16 fm -3 or 2.7 x 10 14 g/cm 3 . At such a high density, nucleons in nuclear matter are crushed into their constituents, 

1. e. quarks and gluons. Itoh proposed a model for quark stars in 1970 Q. The deconfinement transition to quark 
£Sj ■ matter was suggested by Collins and Perry in 1975 J2j based on the asymptotic freedom in quantum chromodynamics. 

| In the same paper they also mentioned the possibility that quark matter could be a superfluid or a superconductor 
resulting from the attractive interquark force in some channels. Barrois, Bailin and Love developed this novel idea and 
studied the unusual variant of superconductivity in quark matter, which we now call color superconductivity HQ]- 
They did their calculations in the framework of weak coupling approach and did not take into account the dynamic 
screening of magnetic gluons which are dominant agents in pairing quarks. Therefore the gap or equivalently the 
transition temperature they obtained are too small to be of relevance to any sizable observables. About fifteen years 
later the color superconductivity had been re-discovered by several groups who found the gap could be large enough 
to bring some real effects in compact stars [5J, |g, [D, la, III- llfj| • [For recent reviews on color superconductivity, see, for 

example, [UlilElElElEllBElEi] 

At asymptotically high baryon density or chemical potential, the color-flavor locked (CFL) phase is favored The 



o 
in 
o 



43 

Oh 
i 



- 1—1 

X 



rigorous weak coupling approach of QCD [a, 9. 20L l2ll l22t l23t VlA l2fll l26l| is justified. [For recent reviews on the QCD 
weak coupling approach in normal dense quark matter or in color superconductivity, see fifiL fiH. E?! 12^. l2flL EfflL l3ll | .] 

The situation is more complicated at lower densities realized in compact stars. In this case, the strange quark mass 



m s is in the range from its current mass ~100 MeV to its constituent mass ^500 MeV. The main effect of m s is 
to cause a mismatch in chemical potentials between strange quarks and light quarks. The CFL phase is therefore 
broken down to less symmetric phases. Moreover the equilibrium and electric and color charge neutrality also 
give rise to a mismatch between the Fermi momenta of the quarks that form Cooper pairs 

If the mismatch in chemical potentials for pairing quarks is large enough, the conventional BCS pairing becomes 
questionable [2a, L^a, Lill lill LiiJ Ui3| and the true grou nd state of quark matter in compact stars has to include 
different, unconventional superconducting states [44. I45L 146L I47L l48l | . 

Single-flavor Cooper pairing is the simplest option for neutral quark matter. Contrary to other unconventional 
pairing mechanisms, it is allowed for arbitrarily large mismatches between the Fermi momenta of different quark 
flavors. Single-flavor pairing in the color anti-triplet channel is possible only in the symmetric spin-one channel 
H I22I I49L IHfl I5IL I52I I53L Im. IHfij . because Pauli principle requires that the wave function of the Cooper pair be 
antisymmetric with respect to the exchange of the pairing quarks. This is in contrast to pairing of different flavors 
where the antisymmetric spin-zero channel is allowed. 

Similar to superfluid Helium-3, where condensation occurs in spin and angular momentum triplets [56], the order 
parameter in a spin-one color superconductor is a complex 3x3 matrix. Therefore various possible phases emerge 
corresponding to different patterns of the matrix order parameter. Among others, there are four main spin-one 
color-superconducting phases: the color-spin locked (CSL), planar, polar, and A phases [HIHHHl- Except for the 



*On leave from Department of Physics, North China Electric Power University, Baoding 071003, P. R. China 



2 



CSL phase, the gap functions of spin-one phases are normally anisotropic in momentum space. The gap in the polar 
phase vanishes at the south and north poles of the Fermi sphere, whereas the gap in the planar phase is anisotropic 
but nonzero in any direction of the quasi-particle momentum. The A phase is special in the sense that it has two 
gapped quasi-particle modes with different angular structures, one of which has two point nodes. Spin-one color 
superconducting phases play important roles in thermal properties in quark matter. Spin-one phases have gaps of 
a few tens of KeVs, which can fit into the cooling data of compact stars The anisotropies especially the nodes 
in some spin-one phases significantly affect various thermodynamical and transport properties like the specific heat, 
the neutrino emissivity, the viscosity, the heat and electrical conductivity, etc.. In Ref. [13], we have computed the 
neutrino emissivity from direct Urea (DU) processes and the specific heat for the mentioned spin-one phases to obtain 
the resulting cooling rates for compact stars. We have discussed the reason why the distribution of neutrino emission 
from the A phase breaks reflection symmetry in space y^|. Similar calculations of neutrino emissivity in DU processes 
in the 2SC phase can be found in Ref. |6(ij . 

We know that DU processes in normal and color superconducting quark matter is controlled by the phase space 
opened up through Fermi liquid properties jfijj, which arises from the quark-quark forward scattering amplitude. 
Up to now the amplitude has been derived only via one gluon exchange (^3, EH , which is a standard weak coupling 
perturbative approach with zero quark mass. In this paper we derive the quark-quark scattering amplitude via one 
gluon exchange but with massive quarks. We also do the same in the NJL model. The forward scattering amplitude 
gives a general formula pf = A*(l — K ) f° r the relation between the Fermi momentum pf and the chemical potential 
/i, where K is the Fermi momentum reduction coefficient with respect to [i and depends on coupling constants. The 
phase space for DU processes with massive quarks both in QCD and in the NJL model is then openned up resulting 
from the Fermi momentum reduction. We find an interesting behavior of k(/j,) as a function of /i in the weak coupling 
regime that «(//) increases with increasing fi and approaches its limit value as fi — ► oo. It means that the phase space 
for neutrino emissions is larger at high baryon densities than at low ones. The property seems robust and independent 
of specific models in computing quark-quark scattering amplitudes and Landau coefficients. 

In this paper we also extend our established formalism 0, -">9| by including quark masses to calculate neutrino 
emissivities in a spin-one color superconductor. The quark masses bring much complexity to quasi-particle dispersion 
relations in spin-one phases, as was shown in Ref. [54l64j- We use the effective quark mass on the Fermi surface in 
the calculation, which is determined by the quark-quark interaction in the Fermi liquid. By assuming a natural form 
of the order parameters for spin-one phases, we show that the quasi-particle dispersion relations bear the same form as 
in the massless case, which much facilitates the calculation. In the collinear and small mass limit, we derive neutrino 
emissivities with massive quarks and identify mass corrections. The mass corrections are of the same order as in the 
massless case. We find no mass correction in the polar phase. The mass correction for the A phase approaches zero at 
asymptotically low temperatures because all contributions are from gapped modes. The mass corrections in the CSL 
and Planar phases are found to be about 26% and 15% of the contributions in the massless case at low temperatures 
respectively. This will have sizable effects on the cooling of compact stars. These percentages for the mass corrections 
are universal as long as n <C 1, which is fullfilled in QCD weak coupling approach and for realistic values of coupling 
constants in the NJL model. Note that mass corrections could be very different near the chiral phase transition or for 
strange quarks with even larger quark masses. However our current formulation is only valid at small mass limit with 
respect to the chemical potential. For very large quark masses, the Landau Fermi liquid theory has to be reformulated 
and the integral in the emissivity is very different from the small mass limit. This is beyond the scope of the current 
work. A future investigation of effects with large quark masses is needed. 

Our convention for the metric tensor is g' 11 ' — diag(l, —1,-1,-1). Our units are h = c = fc^ = 1. Four- vectors are 
denoted by capital letters, K = = (ko, k), while k — |k|. 



II. PHASE SPACE FOR DIRECT URCA PROCESSES IN QUARK MATTER AND EFFECTIVE 

QUARK MASS 

We study DU processes for light quarks u + e~ — > d + v e and d — > u + e~ + V e , where e~ denotes electrons and 
v e and V e denote electron neutrinos and anti-neutrinos respectively. The phase space is essential for these processes 
to proceed. In normal quark matter, if quarks are treated as free and massless, (3 equilibrium requires [id = /tt u + /x e , 
where [id, [in and [i e are chemical potentials for d, u quarks and electrons respectively. This leads to the relation 
in Fermi momenta pFd — PFu + PFe- One sees that the phase space is zero because there is no triangle inequality 
among Fermi momenta. If the quark-quark interaction is switched on, Fermi momenta are not the same as chemical 
potentials any more, instead they get negative corrections from Landau Fermi liquid property (6.2, 63], 
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Figure 1: Phase space for DU processes with massless quarks 
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where as is the strong coupling constant and C'f = {N% ~ -O/C^Vc) m the eigenvalue of the Casimir operator of 
SU C (N C ) group in the fundamental representation with the number of colors N c = 3. The corrections to Fermi 
momenta open up the phase space characterized by the triangle inequality pFd < PFu + PFe- As illustrated in Fig. 
two dashed circles denote the Fermi surfaces for free d and u quarks. They shrink to two smaller solid circles after the 
interaction is turned on. The amount of reduction in Fermi momentum for d quarks is larger than that for u quarks, 
then there is a triangle among Fermi momenta implying non- vanishing phase space. The corrections proportional to 
a s is from the quark-quark forward scattering via one gluon exchange with zero quark mass. 

In this section we will re-analyze the phase space for DU processes by deriving pf as a function of fi in QCD with 
quark masses and some other features not considered in Ref. We also carry out the same task in the NJL 

model which has not been done before. 

The relation between the Fermi momentum and the chemical potential in the Landau Fermi liquid theory can be 
derived from 
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where d g = 2NfN c is the quark degeneracy factor and Nf the number of flavors, n is the number density of quasi- 
particles. The Landau coefficients / j are defined by 
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Here f pa ,p'<r' is the Landau Fermi interaction for quarks with spin indices a and a'. The angle 6 is between p and p' 
both taken on the Fermi surface. The integral is over all directions of p. The factor 1/4 arises from the average over a 
and a'. The factor (21 + 1) comes from the normalization condition for Legendre polynomials J dxP^(x) — 2/ (21 + 1). 
The density of states on the Fermi surface is given by 



N(0) ee I drS(E pa - M ) = M 
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where (dp/ dE v<J ) PF is the inverse of the quasi-particle velocity on the Fermi surface. Note that we have assumed that 
the flavor dependence of the quark mass is weak and negligible, which is true for j3 equilibrium where fx u w fid, so we 
suppressed the flavor indices of the chemical potential fi. 

In QCD, the quark-quark scattering amplitude via one-gluon exchange is given by 
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An average over colors, flavors and spins in the initial state and a sum over all of them in the final state have been 
taken. Here P\ = (-E Pl ,pi) and Pi = (E P2 ,p2) are on-shell 4-momenta of two colliding quarks respectively, m is the 
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quark mass, and g is the coupling constant in QCD and as — g 2 /(47r). The amplitude diverges at small scattering 
angles 9 since (Pi — Pi) 2 w — 2p F (l — cos 6) — ► if m ^ 0. The Landau coefficients can be obtained by applying Eq. 
©• Of course all Landau coefficients are divergent. However the combination f§ — |/f is finite, 



Here the momenta and energies are set to the Fermi momentum and the chemical potential respectively. Substituting 
the above expression into Eq. (J2J, one obtains 
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where we have used n = j^Pp- We assume as is small. The Fermi momentum is related to the chemical potential 
in the form 

p F = p(l- K ), (7) 

where k is the Fermi momentum reduction coefficient with respect to p and is assumed to be a small dimensionless 
number. Then the solution of Eq. (jSJ reads 
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where we only kept the leading contribution. One can see in Eq. © the quark mass term belongs to the next-to- 
leading order and does not appear. The Fermi momentum reduction coefficient n(p) in Eq. {jJJ is illustrated in Fig. 
|2j There are two unconnected branches in k(/x), one is above the static solution k = c ^° s , the other below it. For 
the upper branch, n(p) monotonously decreases and approaches the limit vaule k — c ^ as with increasing p. For the 
lower branch, the trend is opposite. Both branches go towards the limit k — at asymptotically large p, i.e. 

k{/i — > 00) = G o~ s ■ But at lower p or lower baryon densities, there are two different trends for k(h), one corresponds 
to more phase space and the other to less phase space for neutrino emissions. 

If k is small and its dependence on p is weak, Eq. jSJ can be simplified to an algebraic equation without assuming 
as is small. The solution is 

C F a s (g) 



2tt + 5C F a s 



The above solution is lower than the limit value n(p — > 00) = ^f* s and can approach it when as < 1. It is interesting 
to compare the solution © with @. If near a specific value fio the dependence of k on the chemical potential is 
negligible, k(jj,q) is given by Eq. JgJ satisfying k(ho) < c ^° s , which leads to n(p) < f° r an y values of p following 

Eq. 0. Thus the lower branch of the solutions in Fig. [2 is physical, implying quenched phase space for neutrnio 
emissions at lower baryon densities. 

One can also determine from Eq. Q and © the effective quark mass m ~ py2~K which is proportional to the 
chemical potential. Note that the effective quark mass is determined by the interquark interaction on the Fermi 
surface and can be much larger than the current quark mass. 

The above results are from the QCD weak coupling approach and valid at very high densities. At mediate densities 
realized in neutron stars, one has to use more phenonmenological model like the NJL one. For a review of the NJL 
model in quark matter, see, e.g., Ref. ^J- Now we consider the simplest case in the NJL model, a two flavor system, 
whose Lagrangian is 

C = i> [vfdp - m] i/j + G s [(^) 2 + (ferV) 2 ] - G v (^T a ^) 2 , (10) 

where T a = X a /2 and A a are the Gell-mann matrices in color space, r = (t^t^t 3 ) are Pauli matrices in flavor 
space. Gs is the coupling constant for the scalar and pseudoscalar channel, and GV is that for the vector channel. In 
some sense, these two channels are independent of each other. Hereafter we study each channel separately by setting 
the coupling constant zero in the other channel. 

To obtain Landau coefficients /jf ls we calculate the quark-quark scattering amplitude, 
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Figure 2: The Fermi momentum reduction coefficient k(/i) as given in Eq. 0. The dotted line is the static solution k = c ^" s , 
The solid line is the branch of /t(/i) below the dotted line corresponding to an initial value smaller than , while the dashed 

2ff£ ± 0.08; (b) with «(|Lto = 0.3) = 



line is that above it. (a) With the initial value K.(fi = 0.3) 
(JU, the lower branch (solid line) is the physical solution. 
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where B\ and i?2 are two coefficients depending on coupling constants Gs and Gy, 

B x = 2G s + C F Gv, 

B 2 = G s {2N c Nf + l) + 2C F Gy 

The amplitudes As, Aps and Av for scalar, pseudoscalar and vector couplings respectively are given by 



(12) 



A S 
A PS 
A v 



1 



2N f N c 
3 

2N f N c 



G s 



Pi-P 2 - (4N c N f - l)m 2 



E E 



G s 



Pi ■ P 2 - m 2 
E E 



NfN c 



G\ 



Pi-P 2 



2m 2 



(13) 



An average over colors, flavors and spins in the initial state and a sum over all of them in the final state have been 
taken. The Landau coefficients can be obtained by using Eq. JS}, 
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where the momenta and energies are set to the Fermi momentum and the chemical potential respectively. Substituting 
Jq X of Eq. 11411 into Eq. (2J, using n — -^pp, one obtains 
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If we make the assumption Bi t2 fi 2 /ir 2 ~ Bi^ 2 p 2 F /-K 2 ~ k< 1 (actually it is a good approximation for realistic values 
of Bi, but not so good for B 2 ) together with all other assumptions in solving Eq. ©, we get an analytical solution 
for k, 
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which is in a similar form to Eq. JHJ- The solutions of k(^) also have two branches in analogy to the solutions of 

Eq. ||HJ shown in Fig. (2J. We will argue that the lower branch k(^) < is favored, the same situation as in the 

QCD case. 
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If k is small and independent of one can solve Eq. II 1511 without the assumption -Bi^A* 2 / 71 " 2 ^ 1) similar to the 
QCD case, 

(17) 



3tt 2 + (7B 1 + 3B 2 )^ 2 



Since we have assumed that K is independent of fx, B\^^ 2 should be constants. The above solution is smaller than 

the limit value 2 -§^r- of in Eq. (1 1 fil) and recovers it when -Bi^/i 2 / 71 " 2 ^ 1- Therefore we have < 2 
for any values of \i following Eq. iflfijl . Therefore both the QCD and NJL cases give the consistent trend of as 
monotonously increasing functions of towards their limit values at high densities. 

We will use Eq. \\7\ in our forthcoming calculations. The effective mass on the Fermi surface is given bym~ fj,y2n. 
In the scalar and pseudoscalar channel, we set Gg = 5.1 x 10 _6 MeV~ 2 , jti = 500 MeV (these parameters are the same 
as in Ref. [65]), then we obtain k w 0.052 and m w 160 MeV. In the vector channel for Gy — 4.5 x 10~ 5 MeV -2 , one 
obtains k ps 0.13 and m w 250 MeV. In both channels, the effective quark masses are much larger than the current 
mass but much smaller than [i. 

From Eq. JgJ or itTTjl , one also obtains the angles between u quark and electron momenta and between those of u 
and d quarks on the Fermi surface, 

n 2 

cosd ue »s 1 — k, cosd u d w 1 — k - — . (18) 

PFuPFd 

These formula are useful in deriving neutrino emissivities in DU processes. 

A few remarks are in order. All the results in this section apply to DU processes in normal and color superconducting 
quark matter. When calculating the quark-quark scattering amplitude, we do not take into account the fact that 
the chemical potentials for u and d quarks are different, the same approximation is also made in deriving Eq. Q in 
QCD In deriving Eq. (JJ, © and ifTTIl . we have made the approximation that pf ~ M- This is consistent 

to the assumption that the effective quark mass is much smaller than the chemical potential. In this paper we take 
this assumption and use Eq. H17I1 in forthcoming calculations of neutrino emissivities. We find in the weak coupling 
regime the solutions for the Fermi momentum reduction coefficient K as monotonously increasing functions of chemical 
potential in the QCD and NJL cases. Such a trend in chemical potential implies that the phase space for neutrino 
emissions is quenched at lower baryon densities. This property seems robust and independent of specific models used 
in computing the Landau coefficients. In general case where pp ~ [i does not hold, Eq. (JJ is still valid but k has 
to be solved self-consistently. All the results in this section are for light quarks. If one takes into account strange 
quarks, the situation becomes much more complicated and one has to reformulate the whole Fermi liquid theory. This 
amounts to dealing with not only scatterings among strange quarks, but also those between light and strange quarks. 
A new Fermi liquid theory with two different Fermi surfaces is necessary. 



III. QUARK PROPAGATORS IN A SPIN-ONE COLOR SUPERCONDUCTOR 

To obtain the neutrino emissivity in a spin-one color superconductor, we calculate the imaginary part of the W 
boson polarization tensor from the superconducting quark loop, which involves the quark propagator S (K). To obtain 
S(K), we start from S ,_1 (i4r) as follows 

S-\K) = ( % g ) = ( P ^tW°- m p t /70 ^ 70 ) , (19) 
V s 2i s 22 J V 4>fM k - m/7o -mj' 

where Mk — A.^^Jk) J\ and 7oA^j c 7o = A* ; 7^(k)J; are 12 x 12 gap matrices in color and Dirac space, ptf is the 
chemical potential for quarks of the flavor /. 0/ is the magnitude of the gap. The form of An defines a specific 
spin-one phase, see Table I of Ref. [58] for An in various phases. Here 7j_(k) = 7 — k(k • 7) are transverse Dirac 
matrices perpendicular to a momentum direction k. and J/ are generators of SO(3) algebra with elements (Ji)ij — it 11 -' 
with the anti-symmetric tensor As in Ref. y^, we only study transverse phases. 

The structure of Mu in Eq. iflfljl looks like in the massless case, but it is different actually. In the massless case 
S21 is expanded in terms of massless energy projectors A£. m=0 [put m = in Eq. IIAlj) ] as S^ 1 = A£ m _ Q <jffM.\ s _. 
The expansion is justified by the fact that A£ m=Q is commutable with 75 and 7j_(k). In the massive case the energy 
projectors A£ defined in Eq. (|Alj) are not commutable with 75 and 7±(k). So if one still assumes the same form 
of the gap matrix, the propagators and the dispersion relations would be much involved as are manifested in Ref. 

mull. 
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The propagator S(K) is found by taking the inverse of S 1 (K) in Eq. see Appendix El for the derivation. 

One then obtains a very transparent form of the dispersion relations, 

e e KrJ = ^{^f-eEkY + X^), (20) 

where Ek = V k 2 + m 2 . Ak,r denote eigenvalues of the matrix 7oAlj c A / Jk7o given in Table I of Ref. :.58] for various 
phases. The gaps in a specific phase are given by \J\^ r 4>f- This form is the same as in the massless case except 
that k is replaced by Ek- The results for Sn and S22 are given in Eq. l|A2l) and \ A3I) . Since only quasi-particles are 
relevant in neutrino processes, we only keep excitations of positive energy, the branch with e = +, in Sn and S22. 
Then we have 

^11/22 = — T2 ~2 'k.r A k70, 

e k,r,/ 

where we have suppressed the superscript e = + in the quasi-particle energy, £k.r,/ = e k r /- The projectors are 
given m Eq. Q . 



IV. TIME DERIVATIVE OF THE NEUTRINO DISTRIBUTION FUNCTION 

Our starting point to compute the neutrino emissivity is the time derivative of the neutrino distribution function 
derived from the Kadanoff-Baym equation for DU processes u + e~ — > d + f and rf — * u + e~ + 77. The procedure can 
be found in subsection A of section II in Ref. Ififll. The derivation of Kadanoff-Baym equation can also be found in 
various papers using closed time path formalism |6dl6^.lfi3 |. The result is 

d G 2 f d?p 

Ot 4 J (^TTJ PuPe 

Here n F / B (x) = l/(e x ± 1) are Fermi or Bose distribution functions. fv(t,Pv) is the neutrino (including the anti- 
neutrino contribution) distribution function. p e , p„ are electron and neutrino momenta respectively. fi e is the chemical 
potential for electrons. Gp is the Fermi coupling constant. The leptonic tensor L\ a (p e ,p v ) is given by 

Lx,(Pe,P,) = Tr [ 7A (1 - 7 5 ) 7M P ( r 7CT (l - 7 5 )7p^1 

= 8 (ie XaPePu + P£P° + P»Pl - g Xa P e ■ P v ) , (22) 

where P e and P v are on-shell 4-momenta. We have used the notation exap a p v = exa^rP^P^ with e\ avT the anti- 
symmetric tensor. U^(Q) is the retarded self-energy for W bosons with a superconducting quark loop, where 

Q = (Pe - - Pu,Pe ~ Pv)- 



V. IMAGINARY PART OF OF THE W-BOSON POLARIZATION TENSOR 



Having quark propagators, we can calculate the imaginary part of the W-boson polarization tensor. Following the 
same procedure as in Ref. [53, one arrives at the following expression, 



Imn ACT (01 - -Im W 1 U k o-Vu + E k }[p ~n d + E p } Ag 
" \ \ 4E k E p { (k 2 - el r J(p 2 - el s J ,_t 

[kg + - E k ] bo + Pu~ E p ] Xa - 1 
(fco 2 -<^)^o-4. S ,J "'^ J ' 
Here we used P = K + Q. The quark tensor T^f± involves color and Dirac traces 

T r x ° ± (k,p) = 4£; fc i? p Tr[7 A (lT75)P k ± r A±7V(lT75n ± s A±7 



(k,p) 



(23) 



(24) 



The negative component T x f_ can be related to the positive one T x s a + through Eq. Ip8)l. We can rewrite T r x s a + (k, p) 



as 



7# + (k,p) = Tr 7 A (1- 75)^(7^)7^1-75)^(7,^) 



(25) 
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where we have dropped all mass terms which are vanishing. One sees that "f^K^ and ^ V P V appear in T x °, (k, p), 
where K 11 and P v are on-shell quark momenta. In normal quark matter, one of the projectors, say V£ r=1 can be set 
to 1 and others with r ^ 1 set to zero, then one reproduces the result in the normal phase, 

T X °(K,P) ee Tr[ 7 A (l- 75 )( 7M ^) 7 CT (l- 75 )( 7 ^)] 

= 8 (ie XaKP + P x K a + P a K x - g Xa K ■ P) . (26) 

Note that the anomalous propagators S21 and S12 containing the uu and dd condensates do not contribute to the 
self-energy. This is in contrast to spin-zero color superconductors, such as the 2SC and CFL phases, in which the 
anomalous propagators are off-diagonal in flavor space. This difference is related to the conservation of electric charge. 
Performing the Mastubara sum, one can extract the imaginary part, 



n-\a /r -\ne, ne, n F{~ eiew,u)nF (e2£ p ,s,d) r/ > 

T x ° + (k, p)B£ B% d \ p / ' \ S (lo - eie k , r , u + e 2 e p , s , d ) 
+T r _(k,p)B k 1 r d B 2 ^—5{q a + eie ktr>d - e 2 e p , s ,„) 



(27) 



where qa = p e — He — Vv ■ The Bogoliubov coefficients are defined by 

£k,r,/ + e(/i/ - Ek) 



B Kr,f 



with / = u,d and e = ±. The two terms inside the square brackets are charge-conjugate of each other and give the 
same result. To see this, one changes the summation indices e\ «-> 62,?" <-* s in the second term, introduces the new 
integral variable k — > — k q = p (also means p = k + q — > — k), and uses 

r/ r CT _(-p,-k) = r r ^ + (k,p). (28) 

The proof of the above identity for all phases considered in this paper is given in Appendix [5] Taking into account 
Ak.r = A_k,r and then ek,r./ = £-k,r,/, we end up with the first term in Eq. 112 7|l . So we can keep the first term in 
Eq. ll2?j) and double the result, 

/J 3 k 1 „ 



(2tt) 3 4^ fe Sp 
uf{— eiek,r,ii)nF(e2e P ,s,d) 



«s(-ei£k,r,u + e 2 e P , s ,dj 



(5(<?o - eiek, r ,« + e 2 e P)S)d ). (29) 



VI. QUARK MASS CORRECTIONS TO NEUTRINO EMISSIVITY 

In this section we will calculate neutrino emissivities with massive quarks and identify the mass corrections. We 
will use the effective quark mass m ~ /j,-^2k given in Sect. |n| The effective quark mass is formed on the Fermi surface 
by the interquark interaction in the Landau Fermi liquid. Normally the effective quark mass is much larger than the 
current mass, but it is still small compared to the large chemical potential. 

Inserting Imn A i cr ((3) in Eq. ll2*9"|l back into Eq. i|2~Dl . one obtains 

xn F (p e - He)nF(~eie^ r ,u)nF(e2e P ,s,d)S(qo - ejek.r.ti + e2e P)S ,<j), (30) 

where we have replaced the integral over p e with that over p. In normal quark matter, the quark tensor T Xa in Eq. 
l|26l) contracts with the leptonic tensor L\ a in Eq. ll22l) giving the matrix element of the DU processes 



T x °L XrJ = \M\ 2 , 



(31) 
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with \M\ 2 = 256(P e ■ K)(P V ■ P). Using Eq. (JJ and JEJ, the matrix element |M| 2 can be worked out, 

|M| 2 « 512^^^ [1 - (I-k)cos^], (32) 

with k given in Eq. ifTTjl . 

The quark tensor ^^(k, p) in spin-one phases are evaluated in AppendixJO Using Eq. IIClll . its contraction with 
L XcT is given by 

L Xa T Xa + = u; rs \M\ 2 + L Xa 5T x f . (33) 

Here the values of u rs are the same as in the massless case (see Tab. II of Ref. [^| or Appendix EJ. As is 
shown in AppendixEl the dominant contribution in 8T x s a is from the term proportional to RkT X(T (Ko, P) with Rk = 
(Ek — k)/k m 2 / (2p 2 Fu ) ~ n and Kq = [k, — k), because P e -Ko ~ 2p e fc is much larger than P e -K ~ 2nE k p e . So f>7^ CT 
can be written as (57^ CT w Xr S RkT Xo (K 0i P) where Xrs is listed in Tab. □ Terms of R p T Xa (K, P ) ~ nm 2 /p 2 Fd ~ k 2 
are next to leading order since the appearance of Po does not enhance the order of magnitude of the result, so they 
are suppressed by k compared to the massless term. Also terms of R k R P T Xa (Ko, Pq) ~ to 4 / {p 2 Fu p 2 Fd ) ~ k 2 are also 
suppressed due to double R factors. Explicitly one derives 

2 

L Xa 5T Xa « Xr S — 256E k E p p ePu {l-(l- K )cos9 du }, 

~ XrsT^—256E k E p p e p u [1 - (1 - k) cos6» d „] . (34) 

The term ncos9 dv inside the square brackets can be dropped since it is of high order. 

Let us understand Tab. [I] We find Xrs = for the polar phase. The reason is simple: all projectors ilC2j) are in color 
space or decoupled from Dirac space, so the quark tensor T r Xa + is proportional to T Xa . For other phases, A, planar 
and CSL, color and Dirac indices are coupled. Thus with massive quarks, the Dirac structure of the quark tensor is 
more complicated. Taking the A phase for example, it is easy to understand Xrs — for rs = 13, 31, 23, 32, 33, because 
the excitation branch 3 only involves the color blue while branch 1 and 2 involve colors red and green. One observes 
that X12 = X21 7^ even in the collinear limit, meaning that there is an entanglement between two different branches. 
Remembering that the projectors of the A phase can be expressed in terms of helicity projectors, as shown in Eq. 
(60) of Ref. quasiparticles of branch 1 have helicity +1 for k 3 < and helicity -1 for k 3 > 0, while quasiparticles 
of branch 2 have opposite helicities. In the massless case, helicity means chirality. The DU processes only allow 
quarks with left chirality to participate. Therefore in the collinear limit fc 3 w p 3 , different branches cannot crosstalk. 
However when quarks have masses, a helicity eigenstate is not chirality one. In an eigenstate of left chirality, one can 
find eigenstates of both helicities, which makes two different branches coupled together. In the same way, one can 
understand the appearance ST r Xa in the quark tensor with different structure from T Xa . A projector of a spin-one 
phase can be expanded with respect to helicity projectors. In the massless case, the quark tensors for all branches 
have the same Dirac structure as in normal quark matter. With non-zero quark mass, this is no longer the case: 
new terms in the quark tensor emerge from an overlapping of subspace with different helicities (but with the same 
chirality). They have different structure and cannot be collected into T Xa as in the massless case. This can be seen 
from the fact that such terms as 707 • h'y^K^ or j^K^^f ■ k7o appear inside the trace of Eq. Il2.'it for A, planar and 
CSL phases, 

707 • £7^ = 7^(7 ' k) 7 o = 1^ - Rk7„Kg. (35) 

When masses are zero, the result is just 7^^^, which results in traces proportional to T Xa in Eq. il26|) . When masses 
are non-zero, there is an additional term Rkj^K^ different from T Xa . 
Using Eq. 13011 and 13411 . the emissivity is written as a sum of two terms 

= - J ^„|/,(t,p,) = ^+eW. (36) 

Here is in the same form as the emissivity in the massless case which we denote as ei™ -0 "* , i.e. Eq. (40) of Ref. 

except that all k and p are replaced by E k and E p in dispersion relations respectively, and that 2as/(37r) replaced 
by K in Eq. © or itTTjl . It comes from the term w rs |M| 2 in Eq. In the massless case, the ranges of the integrals 

over fc — n u and/or p — jid on which the quark distributions depend are [— fi u /d, 00], which are further taken to be 
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Table I: The values of Xrs in quark mass term 5T r x / « Xr S RkT Xa (Kq, P)- For the polar phase Xrs = for all r and s. The row 
and column are labeled by r and s respectively. The collinear limit is taken. We used the notation £ = fc 3 and 9(0 = 1 +sgn(£) 
for the step function of £. 



Planar 


1 


2 


1 


2£ 2 


2£ a 






2 


2£ 2 


2£^ 


d+« 2 ) 2 


d+e 2 ) 2 



A 


1 


2 


3 


1 


-o(0 







2 


0(0 


-0(-O 





3 












CSL 


1 


2 


1 


i 

2 


i 
2 


2 


1 
2 


1 
2 



[—oo, oo]. With non-zero quark mass, the integrals over k — fj, u and/or p — [id have to be transformed to those over 
Ek/p — Hu/d on which the quark distribution functions depend, i.e. 

dpdk = (E p E k /pk)dE p dE k = (E k E p /pk)d(E p - fi)d(E k - jj) 



pk 



-d(E p - n)d{E k - n). (37) 



In the last equality we have chosen E k « \i u and E p w fid on the Fermi surface. The range of E k / p — \i u j d becomes 
[m — fi u /d, oo] and is further approximated to be [— oo, oo] if m <C fJ. u /d- Using Eq. itlSll . (13771) and l(33j) . we obtain 

e (o) = e (m=o) (38) 

This means ei is identical to the emissivity in the massless case. 

The second part of the emissivity in Eq. ijflfif . e£ , has obvious mass dependence. It comes from the term L\ G 5T^f 
in Eq. l(33j) . Inserting Eq. ll^ljl into Eq. l|30fl and using Eq. pfif. one obtains 



k 2 /" d 3 p„ f dE p dn p f dE k dn k v-^ , . 

32 -— G F/Wd 7 (2^ 7 "W 7 l^F^ ^ X-(l-cos^) 



(2tt) 3 ^7 (2tt) 3 J (2tt) 3 

xB k,r,u B p,s,d n F(Pe ~ Ve)n F (-eiek, r ,u)n F (e 2 £p, s ,d)£ [COS 6»„rf - 1 + K[ll/(pFuPFd)\, (39) 

where Eq. lj3Tjl was used in the integral. We have also used 

5(p e - fi e ~p v - eie kj . r!U + e 2 e p . s . d ) « — — — <5[ cos _ i + Kfi 2 e / (p Fu p Fd )}. (40) 

PFuPFd 

After taking the same procedure as in Ref. j^|, one has 

4 X) « —57— — GF^e^u^dT 6 G' (ip u , (fd), (41) 

with integral G' defined by 

/l />oo />oo />oo , 

d£,Xrs(0 dvv 3 dx dy n F (v + e^y 2 + X^Ofl - ^Jx 2 + X s (0fd) 
-1 Jo Jo Jo 



ei,62 r,s 



xn F (- ei ^/y 2 + \ r (Z)<pl)n F (e 2 y/x* + A s (£)^)- (42) 

Here we used dimensionless gaps scaled by the temperature, ip u _d = 4>u,dlT . Eq. (14 XII shows that el. 1 '' is in the same 

order as ei m_ °\ We can take the special case ip u — ipd — f and do the numerical calculation for G' as a function 
of ip. The results are shown in Fig. <PJ. We see that G < for all three spin-one phases. At asymptotically low 
temperature or large <p, G' for CSL and Planar phases approach constants due to the gapless mode 2, 

G' = (CSL), - 45? - 6 ^ 2 - 1 ) (Planar). 

5040 v h 5040 v ' 

For the A phase, G' turns to zero at large ip since all contributions come from the gapped modes. 
In summary, from Eq. ffify. (l38)l and I^U, the emissivity can be written as the sum of two terms, 

e I/ = e ™=° + 4 1) . (43) 
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where e„ is the correction due to non-zero quark mass and it is of the same order as t v ■ It is zero for the polar 
phase. For the CSL and Planar phases, we obtain 

4 1} « Ci— ^4 m=0) ~ C^ m =°>, i = A, CSL, Planar 

1 — K 

At asymptotically large ip, the coefficients C'i are given by 

C A = 0, C CS L = -i«-0.25, C Planar = -^-^ m-0.15. (44) 

which show that mass corrections are about 25% (CSL) and 15% (Planar) of the contributions in the massless case at 
very low temperatures. These results are quite universal at k -C 1, which is the case for the weak coupling approach 
in QCD with one gluon exchange as n ~ as- This is also the case for the realistic values of coupling constants 
G s ~ 5.1 x 1CT 6 McV~ 2 and G v ~ 4.5 x 1CT 5 MeV -2 in the NJL model which give k ~ 0.05 - 0.13. 

A few comments are in order. The gap matrix in Eq. II19JI leads to the same form of the dispersion relations as 
in the massless case. This much facilitates the calculation with transparency. If the gap matrix is not in this form, 
the dispersion relation is more involved 0,0]. But one can P rove that they are equivalent in some limits (see the 
discussion in section VII of Ref. We have used the collinear limit, where the momentum direction of u quarks 

is almost the same as that of d quarks. The effective quark mass m ~ is also used in the calculation, which 

is much larger than the current mass but small compared to the chemical potential. Then we find that the mass 
corrections to the emissivity are of the same order as in the massless case. The mass corrections from gapless modes 
in CSL and Planar phases which are dominant at low temperatures are about 25% and 15% of the contributions in 
the massless case respectively. These percentages for the mass corrections are universal as long as k C 1. The mass 
corrections could be very different near the chiral phase transition and for the DU processes with strange quarks with 
even larger quark mass. However the current formulation is only valid at small mass limit compared to the chemical 
potential. For very large quark mass, the Landau Fermi liquid theory has to be reformulated and the integral in the 
emissivity is very different from the massless case since the integral range E — \x cannot literally be taken as [—00, 00] 
due to m ~ //. Not only this has to be modified, one also has to deal with scatterings between light and strange 
quarks, which involves the Fermi liquid theory with two different Fermi surfaces. This is beyond the scope of the 
current work and deserves a future study. 
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Appendix A: ENERGY PROJECTORS AND PROPAGATORS 



We make use of energy projectors to get S(K) from S (K) in Eq. I|19ll . The energy projectors are given by 

7o(7 • k + m) 



1 + e- 



(Al) 



where Ek — \Jk? + m 2 . We can express and S 2 2 m terms of energy projectors: 

S Ti/22 = ( fc o ± M/ - eS fe )7oAk. 
The chemical potential (if has flavor index /. The elements of S(K) are given by 

Sn = [S ± i — S' 12 1 (5 22 1 ) 1 S' 21 1 ] , 

5*22 = [S 2 2 — SziiS^) 1 <5i2 1 ] 7 

5*21 = —(S22) S21 Su, 

S12 = —{S11) 1 S 12 S22- 

To get Sn, we need to evaluate 

Here we have used A k aj_(k) = aY(k)A k e with a = 707. Note that 

L+=[ai(k)ASJj[ai(k)Ai/rJi/] 

is just the matrix 7o-Af k .Mk7o in Eq. (14) of Ref. and it is commutable with A k . We can decompose in 

terms of its projectors as = \k,r'P^ r . Then we obtain 



k - (if - eE k 



A: 



fc - /// + eEfe , 



^r A k70, 



\2 k,r 



Sn 



with the quasi-particle excitation energy e k r ^ = J (fi — eEk) 2 + ^k,r<P 2 f 
Making the change (if — > — /i/ in 5n, we get S22 

fc + (if - eE k , 



(A2) 



>J22 — 



^0 ( £ k,r,/) 



2 _77 k,r A k e 70, 



(A3) 



where P k are projectors of 



L k = [a' ± (k) Ay J/] [ai (k) A *, V J V \, 
whose eigenvalues Ak,r are the same as Lt's. The projectors V^ r are given by 



k./ 



TT L k ~ ^k,r 



(A4) 



Note that the projectors V^ r are commutable with A k since [L k ,A k ] = 0. 



Appendix B: PROOF OF THE IDENTITY (J28j| 

One starts from T^f_ as follows 

7#_(-p, -k) = Tr [ 7 A (1 + 75)^r p , s A: p7 ° 7 CT (l + 75)Pl k , r A: k 7 



Inserting the charge conjugate operator C = ijoj 2 and using CoH'oPC 1 = (aPot l ) T , a\(— k) = ct l j_(k) and J; 
one obtains 

One uses C-f^C^ 1 — to obtain 

CAlp/C- 1 = (A+ 7 °) T 
Also one has C75C" 1 =75 =75'. Then one finally arrives at 

T£_(-p, -k) = Tr [ 7 AT (1 + 75) T (^ p + , J T (A+ 7 °) T 7 ctT (1 + l,f W^Y '(A+ 7 f 

= Tr [ 7 A (1 - 75)A+ 7 X r 7 CT (l " 75)A+ 7 °P+r 
- Tr [ 7 A (1 - 75)^ k +,A+ 7 7^(l - 75)Pp+ s A+ 7 ° 
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^r+(k,p), 



where [P£ r , AJ7 ] = has been used. 



Appendix C: CALCULATE QUARK TENSOR T£°+ FOR SPIN-ONE PHASES 

The quark tensor T$£+ can be written as a sum of the leading term with the same structure as in the massless case 
and the mass correction term, 



Act / 



-Act 



(CI) 



where T A<T (K, P) is the quark tensor in the normal phase and given in Eq. (126)) . ST r x f is the correction from the 
quark mass. 



1. Polar phase 

The polar phase is particularly simple. The order parameter or the condensate has the form 

M k = J 3 7-L,3(k). 
The projectors V£ r do not depend on the quark momentum k, 

' k,l — J 3 ! 



k.2 



i-Ji 



(C2) 



The trace can be decoupled into a color and a Dirac one. The color trace is 



k,2' k,2J 
+ -D + 



0. 



The Dirac trace is the same as in the normal phase in Eq. l|2fijl . So we have 

= 2, U22 = 1, W12 = u>2i = 0, 
<5T r A / = 0(r, s = l,2). 
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The order parameter is in the form 



2. Planar phase 



M k = j l7 i(k) + j 27 i(k) 



The projectors are 



V + - 



k,2 



1 



1 + kj 
1 — P + 

M 



J 2 (l - fc 2 ) + J2(l _ k \) _ { J 1? J^fo + J 3 fc37 07 5 7 • k 



We can rewrite i and V£ 2 m ^ ne f° rm 



k,2 



1 - Ilk - /2k7o7 7 ' k > 



with 



We obtain 



/ik 
I 2k 



1 + fcf 



J 2 (l - fc 2 ) + J|(l - fc 2 ) _ { J lf J 2 } kl k 2 



1 + fcf 



Ct>22 



= Tr c [Jik/ip] + Tr c [hkhp] , 

= -1 + Tr c [Jik lip] + Tr c [hkhp] , 

= cj 2 i = 2 - Tr c [hwhp} - Tr c [/ 2 kl2 P ] , 

= Tr c [J 2k / 2p ] \-R p T x °(K, P ) - R k T x °{K , P) + R p R k T x °(K , P Q 



5% 



Act 
22 



^ii ff 



Act 



-5% 



Act 



where Rk = {Ek — fc)/fc , K = (J5k,k) and -Ko = (A;, — k). T Acr is given in Eq. l|2fil) . In the collinear limit, k || p or 
k w p, the color traces can be simplified, then we obtain 

CJn = 2, W 2 2 = 1) ^12 = ^21 = 0, 



8T Xa 
ST 22 



n 

Act 



(l + fc§) 2 



-R P T X °(K,P ) - R k T^(K ,P)+R p R k T^(K ,P ) 



5% 



Act 



Act 



-^21 



Act 



The order parameter is 



3. A phase 



7 i(k)+ i7 i(k) 



The projectors are 



V+ 

k,l 
V' 



k.2 



k,3 



2 Js 

2 Ja 
1 - J. 



1 + sgn(fc 3 )7o7 5 7 • k 
1 - sgn(fc 3 )7o7 5 7 • k 



3 i 



where sgn(cc) is the sign of x. The results for u> rs are 



The results for STjf are 



wn = ^[1 + sgn(fc 3 )][l + sgn(p 3 )], 



U>22 



^33 



[1 - sgn(fc 3 )][l - sgn(p 3 )], 



Wl3 — W31 — CJ23 — ^32 — 0, 

W12 = ^[1 + sgn(fc 3 )][l -sgn(p 3 )], 
W21 = i[l-sgn(fc 3 )][l + sgn(p 3 )]. 



--sgn(p 3 )[l + sgn(fc 3 )]P P T ACT (if, P ) 



--sgn(fc 3 )[l + sgn(p 3 )]R k T^(K ,P) 



+ ^sgn(p 3 )sgn(k 3 )R p R k T A(7 (K , P ), 



Act 



-sgn(p 3 )[l - sgn(fc 3 )]i? p T Aff (if, P 



+ -sgn(fc 3 )[l - sgn(p 3 )]P fc r ACT (^ ,P) 



5T& 



+ -sgn(p 3 )sgn(k 3 )R p R k T^(K , P ), 
= 5T A / = 57^ = <5T 2 A / = 5r£ = 0, 
= isgn(p 3 )[l + sgn(fc 3 )]P p r ACT (X, Po) 



--sgn(fc 3 )[l - sgn(p 3 )]P fe P ACT (^o,P) 



<5T 2 Act 



--sgn(p3)sgn(fc 3 )P P P fc r ACT (^o, Po), 
-i S gn(p 3 )[l - sgn(fc 3 )]P p r Aff (K, Po) 



--sgn(fc 3 )[l + sgn(p3)]P fc r ACT (^o, P) 



In the collinear limit k w p, we have 



and 



--sgn{p 3 )sgn(k 3 )R p R k T^(K , Po). 



wn = 29(k 3 ), uj 2 2 = 20(-fc 3 ), lu 33 = 1, 

W12 = W21 = Wi 3 = W 3 i = Cl>23 = w 32 = 0, 



6%Y = 



-6(k 3 )R p T^(K,P ) - 6{k 3 )R k T Aa {K ,P) + -R p R k T (K , P ), 



5r 2 f = -#(-fc 3 )p P T Aff (k, p.) - 0(-fc 3 )P fc r ACT (ifo, p) + -R p R k T (k , p ; 



5% 



AfT 



(JTjf = 0(k 3 )R p T^(K, P ) + 0(-fc 3 )P fc T ACT (if o , P) - -RpRkT (Kq, P ), 



5T 2 f = 6{-k 3 )R p T^{K, P ) + 6{k 3 )R k T x °{K , P) - -P P P fe T ACT (if , Po), 



where 9{x) = 1 + sgn(x) is the step function. 



4. CSL phase 



The order parameter is 



Mi 



J-7±(k)- 



The projectors are 



k,2 



It is convenient to rewrite these projectors in terms of color indices, 



(K2U 



We obtain 



ton 
W12 



1 

4 

W21 = 1 - ^ 



l + (k-p) 2 + 2(k-p) 



1 + (k • p) 2 + 2(k • P ) 



l + (k-p) 2 + 2(k-p) 



and 



6% 



Act 
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(k • P ) -R p T Xa (K, P ) - R k T A °(K , P) + R p R k T Aa (K , P 



-Act/ 



5% 



Act 



-6% 



Act 
12 



svt-Act 

ov 21 . 



In the collinear limit, we have 



wn = 2, w 2 2 = 1, W12 = w 2 i = 0, 



Act 



1 



-R p T^(K, P ) - R k T^(K ,P) + R p R k T (K , P ) 
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